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Abstract 

Within the context of a first-order phase transition in the early Universe, we 
study the collision process for vacuum bubbles expanding in a plasma. The 
effects of the plasma are simulated by introducing a damping term in the 
equations of motion for a U{\) global field. We find that Lorentz-contracted 
spherically symmetric domain walls adequately describe the overdamped mo- 
tion of the bubbles in the thin wall approximation, and study the process of 
collision and phase equilibration both numerically and analytically. With an 
analytical model for the phase propagation in 1+1 dimensions, we prove that 
the phase waves generated in the bubble merging are reflected by the walls of 
the true vacuum cavity, giving rise to a long-lived oscillating state that delays 
the phase equilibration. The existence of such a state in the 3+1 dimensional 
model is then confirmed by numerical simulations, and the consequences for 
the formation of vortices in three-bubble collisions are considered. 
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I. INTRODUCTION 



According to the standard model and its extensions, symmetry breaking phase transitions 
are expected to occur in the early universe. The mechanism by which these transitions may 
take place could be either by formation of bubbles of the new phase within the old one 
(i.e., first order phase transition) or by spinodal decomposition (i.e., second order phase 
transition). Although the verdict is still pretty much in the air as to which one would have 
actually taken place, in the particular case of the electroweak phase transition common 
opinion inclines more towards the first of the two possibilities. As this scenario would have 
it, bubbles of the new phase nucleated within the old one (the nucleation process being 
described by instanton methods as far as the WKB approximation remains valid and 
subsequently expanded and collided with each other until occupying basically all of the space 
available at the completion of the transition. The caveat introduced by the word "basically" 
is important though: in the process of bubble collision, the possibility arises that regions 
of the old phase become trapped within the new one, giving birth to topologically stable 
localized energy concentrations known as topological defects (for recent reviews see refs. p| 
) , much in the same way in which these structures are known to appear in condensed matter 
phase transitions. From a theoretical point of view, topological defects will appear whenever 
a symmetry group G is spontaneously broken to a smaller group H such that the resulting 
vacuum manifold M = G/H has a non-trivial topology: cosmic strings for instance, vortices 
in two space dimensions, will form whenever the first homotopy group of M is non-trivial, 
i.e., 7ri(M) 1. 

To see how this could happen in detail, let's consider the Lagrangian 

£ = - 1/($) (1) 

for a complex field $. If is a potential of the type V = f (|$P — "Z]^)^ and its parameters 
are functions of the temperature such that at high temperatures $ = is the only minimum 
of V, while at zero temperature all the |$| = ?7 values of the field correspond to different 
minima, then the structure of the vacuum manifold will be that of S^. As vri(5'^) ^ 1, 
cosmic strings should then be created in this model during the phase transition. 

The way in which this would actually take place is via the Kibble mechanism 0. The 
basic idea behind it is that when two regions in which the phase of the field takes different 
values encounter each other, the phase should interpolate between this two regions following 
a geodesic in the vacuum manifold. In the context of a first order phase transition, a 
possible scenario for vortex formation in two space dimensions would then look like this: 
three bubbles with respective phases of 0, 2tt/3, and 4tt/3 collide simultaneously. According 
to the so-called geodesic rule then, if we walk from the first bubble to the second , then from 
the second to the third, and finally back from the third to the first again, the phase will have 
wounded up by 27r in our path, having traversed the whole of the vacuum manifold once 
along the way. Continuity of the field everywhere inside the region contained by our path 
demands then that the field be zero at some point inside of it, namely the vortex core. In the 
limit in which the bubbles extend to infinity, outwards from the center of collision, removal 
of this vortex would cost us an infinite amount of energy, since it would involve unwinding 
the field configuration over an infinite volume. The vortex is thus said to be topologically 
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stable. In three space dimensions, the resulting object would obviously be a string, rather 
than a vortex. 

Clearly, there are other ways in which strings could be formed. Collisions of more than 
three bubbles could also lead to string formation, or, for instance, two of the bubbles could 
hit each other first, with the third one hitting only at some later time while the phase is still 
equilibrating within the other two. This event in particular will be far more likely than a 
simultaneous three way (or higher order) collision, and it is probably the dominating process 
by which strings are formed (especially if nucleation probabilities are as low as required for 
WKB methods to be valid). If this is the case, the importance of understanding the processes 
that take place when two bubbles collide, as well as the need to clarify the role of the third 
bubble and the conditions that should be met in order to get vortex formation, becomes 
clear. 

Bubble collisions have been studied by a number of groups, most notably by Hawking et 
al. 1^, Hindmarsh et al. Srivastava |^ and more recently by Melfo and Perivolaropoulos 
0, but important aspects of the question remain unclear so far. For instance only in 0] and 
1^ the posterior collision of a third bubble is treated (the first two references studied only two 
bubble collisions) , and in none of the referred works did the authors concern themselves with 
the interaction between the bubble field and the surrounding plasma. From current work in 
the electroweak phase transition however, we can expect that in many cases such interaction 
will not be negligible (e.g, see Ref P). The basic underlying reason is relatively simple to 
understand: as the bubble wall sweeps through an specific point, the Higgs field $ acquires 
an expectation value, and the fields coupled to it acquire a mass. Thus, particles with not 
enough energy to acquire the corresponding mass inside the bubble will bounce off the wall 
(thus imparting negative momentum to it), while the rest will get through. Obviously, the 
faster the wall propagates the stronger this effect will be, since the momentum transfer in 
each collision will be larger, and thus a force proportional to the velocity with which the 
wall sweeps through the plasma appears. 

The aim of this paper will thus be to explore two bubble collisions in an environment 
that damps the bubble motion to try and establish how phase diffusion and equilibration will 
proceed in this case, both analytically and numerically. The damping that we will consider 
here will only affect the modulus of the field though (which is equivalent to damping only 
the bubble wall motion), since, the phase being a Goldstone boson, we do not expect it to 
couple to the plasma as strongly as the modulus. We will find significant differences with 
the undamped case, the main new features being related to the possibility of the phase 
wave catching up with the bubble wall and bouncing off it. Once that is understood, we 
will study collisions with a third bubble and see under what circumstances we can expect 
to have vortex formation. All this will be done in the context of a global U{1) symmetry 
breaking model. 

The paper is organized in three sections: in the first one we establish the general ana- 
lytical model, in the second we study two bubble collision and phase bouncing in 1+1 and 
2+1 dimensions (the 3+1 case should be analogous to the latter one), and finally the third 
is devoted to three bubble collisions in 2+1 dimensions. 
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II. ANALYTICAL MODEL 



Consider the Lagrangian (0) for a complex field $. We will use the same form of potential 
that was used in that is, 



V = \ 



(|<l>| -7]) rj\^\ 



2 VI I 3 



(2) 



This is just a quartic potential with a minimum at |$| =0 (the false vacuum), and a set 
of minima connected by a U{1) transformation (true vacuum) at |$| = ptv = |(3 + e + 

y^(3 + e)^ — 8), towards which the false vacuum will decay via bubble nucleation. It is the 
dimensionless parameter e that is responsible for breaking the degeneracy between the true 
and the false vacua. 

The equations of motion for this system are then 

di,d^^ = -dV/d^. (3) 

For the potential (H), approximate solutions of (|^) exist for small values of e, the so-called 
thin wall regime []TU|], and are of the form 



|$| = ^ 



tanh (x — Rq) 



(4) 



where Rq is the bubble radius at nucleation time and = I ^ P — t^- The bubble then 
grows with increasingly fast speed and its walls quickly reach velocities of order 1. We are 
interested however in investigating a model with damped motion of the walls due to the 
interaction with the surrounding plasma. In order to model this effect, we will insert a 
frictional term for the modulus of the field in the equation of motion, namely 

df.d"^ + 7 |$| e*^ = -dV/d^, (5) 

where |$|= d\^\/dt, 9 is the phase of the field, and 7 stands for the friction coefficient (which 
will as a matter of fact serve as parameter under which we will hide our lack of knowledge 
about the detailed interaction between the wall and the plasma). We have found that there 
are approximate analytical solutions to (^) corresponding to one bubble configurations that 
have reached their terminal velocity in the thin wall limit . Solutions will be found in a 
constructive way: first we will suppose that solutions for which the wall has the form of a 
traveling wave do exist to find out what their terminal velocity would be, and then use this 
expression for the velocity to obtain the detailed analytical form of such solutions in the 
thin wall limit. Comparison with numerical simulations will confirm the existence of such 
solutions. Writing $ in polar form: 

$ = pe^\ (6) 
the equation of motion (^) takes the form 
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d^.d^'p + 7 P = -dV{p)/dp + 2pd^ed^e, (7a) 
{pd^e} = 0. (7b) 

For a single bubble configuration we take ^ to be a constant (the phase of the bubble). 
Equation ( [75|) is then automatically satisfied and the equation left for the modulus of the 
field is 

d^d>'p + ^p=-dv{p)/dp. (8) 

To get an approximate expression for the terminal velocity of a bubble under this equation 
of motion in the thin wall limit, note that because the wall thickness is much smaller than the 
radius of the bubble, we can go to 1 + 1 dimensions. Inserting then the ansatz p = p{x — xo{t)) 
leads to 

(1- xl)p" + (xo +7 Xo)p' = dV{p)/dp, (9) 

where p' = dp/dx. We then get an effective equation for the wall motion simply by multi- 
plying by p' and integrating over — oo<x<+ooto get rid of the extra degrees of freedom 
pertaining to the field, whence 



[Xq +7 Xq) 



/+00 \ r+oo 

p'^dx) = / V'dx = AV, (10) 



where AV is just the potential energy difference between the false and the true vacuum 
phases, and we have used the fact that derivatives go to zero far from the origin. The 
solution to ( pUD for the initial conditions xo{t = 0) = i?o, {t = 0) = is 

1 a 

xoit) = -at+—{e~^'-l) + Ro (11) 

7 7 

where a = AV/ (/ p''^dx). Thus, for values of t ^ 7^^ the bubble walls will have reached 
their terminal velocity 

AV 

= ^(um- '''' 

Although obtained in a slightly different way, this expression coincides with what Heckler 
found in |lT|. Typical estimates for the terminal velocities of thin wall bubbles in the 
electroweak phase transition give Vter ~ 0.1 [§] (although in recent work higher values have 
been found, see 0). 

To get an approximate expression for p valid within this regime , it suffices to rewrite (H) 
with a p = p(r — ro(t)) ansatz, where r is the usual radial coordinate. Using ro= 0, ro~ Vter, 
we get 

(1 - <r)^ + 7--) = dV{p)/dp, (13) 

r being the radial coordinate in 3 spatial dimensions. According to ([l2|) however the terminal 
velocity goes roughly like 
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where 6m is the bubble wall thickness and ptv is the true vacuum value of the field. That 
is, at the values of r for which the first derivative of the field is important (r ~ for a 
thin wall bubble), we have R'jvter ~ (^m <^ 1, and the second sumand in the parenthesis is 
negligible when compared to the first. Furthermore, since the radius of thin- walled bubbles 
is very large, we can also neglect the term (2/r) dp/dr -this is just the standard thin wall 
approximation. Thus, we are left with 



:^-vL)^ = dV{p)/dp, (15) 



whose solutions for the potential written in (0) will be 

{r - Vtert - -Ro) 



Ptv 



tanh 



1 ~ ^ter 



(16) 



which is simply a Lorentz-contracted, moving domain wall. 

We have followed closely the analytical study with numerical simulations of the processes. 
This was done in two steps (for details see Ref 0): first we find a numerical solution of 
the Euclidean equations of motion for the field, whose analytic continuation into Minkowski 
space will give us the shape of the instanton solution immediately after the tunneling has 
taken place. The initial bubbles thus found were then placed in a two dimensional lattice, 
and evolved according to (|^) with a modified leapfrog method and refiecting boundary 
conditions. In the following sections, we will be comparing the analytical results with these 
simulations. 



III. TWO BUBBLE COLLISIONS: PHASE PROPAGATION AND BOUNCING 

We will now try to take a look at the events that take place when two damped bubbles 
collide. To do so, we will organize the section into three parts. The first one will deal with 
the configuration and evolution of the phase from the time when the bubbles are still very 
far away from each other to the point when they have completed the collision process. The 
second one takes on from there to study the propagation inside the bubbles of the phase 
waves that originate at that time, and stretches out until the moment when these waves 
catch up with the walls at the other end of the system. Finally, the third part studies 
the interaction between the phase waves and the bubble walls, finding how the former will 
bounce off the latter and propagate again into the interior of the bubbles, which will thus 
behave almost resonant cavity. 



A. Initial phase configuration 

Within the WKB approximation the bubble nucleation rate is small, and therefore typical 
bubbles are nucleated a long distance apart from each other. The approximate solution for 
a system of two bubbles is then simply the sum of the two independent bubble solutions. 
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^{huhhlel + huhhle2) = pe'^ ~ ^(bubblel) + ^{huhhle2) = pic'^^ + p2e'^\ (17) 

to exponential accuracy. From this it is easy to compute how the phase interpolates between 
the two bubbles 

Pi tan 6*1 + p2 tan 6^2 
tanfe' = , (18) 

which has a kink like shape centered at the origin. One could worry however that, as the 
bubbles move towards each other, phase waves could be generated in the central region that 
would subsequently propagate over the exponentially small (but finite) background of the 
modulus so that (|18|) would no longer adequately represent the combined phase of the two 
bubbles in motion, even if (0) still provided a good ansatz for the combined modulus. To 
check whether this is indeed the case we can solve the equation of motion for the phase (7b) 
using the modulus given by (0). In a 1 + 1 dimensional approximation, near |x| = and far 
away from the centers of each bubble we can take their moduli to be (we will write simply 
V for Vter now) 

Pi ~Pi^e(^+^*-^o)/'''", (19) 
P2 - ptoe-(^-''*+^°)/^'", 

where we have taken the centers of the two bubbles to be initially situated at ±a;o. Then 
(P^) yields for the combined modulus 

p2 ~ p2^e2(''*-^o)/^'- {2 cosh(2x/5„) + 2 cos(^i - ^2)} • (20) 

We can now use this in the equation of motion for the phase ([7E|), to get 

e-e" ^ — e ^ ^ 7 & = 0. 21 

dm cosh(2x/0mj + COS^t^i — U2) 



Note that for an initial phase difference 6i — 62 = n the modulus of the field (20) is zero 
at the midpoint between the bubbles. In this case then the denominator in the last term 
of ( pID goes to zero, which means that the phase has the shape of a step function as we 
go over the origin, switching disco ntinuously from 61 to 62. On the other extreme, if the 
initial phase difference is zero there is of course no dynamics to it. We will then focus in an 
intermediate situation, and find solutions to (0) for a phase difference of 7i/2. Using the 
ansatz 9 = T{t)X{x), we can separate variables to obtain: 

2v 

f+ — T+k^T = 0, (22a) 

Cm 

2 2r 

X" + — tanh(-^)X' + ex = 0, (22b) 

k being the separation constant. The equation for T is obviously that of a damped oscillator. 
For k = 0, a T = const, solution exists consistent with our boundary condition that 6 goes 
to 61^2 as X ^ ±00 at all times, while, for k 0, solutions with temporal dependence on the 
phase will die away on time scales of the order of Sm/v (as a matter of fact these solutions 
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would seem to have been artificially introduced by the switching on of the interaction between 
the bubbles). To exponential accuracy then, the phase will rapidly settle in a stationary 
state that will interpolate between its two values at infinity. The solution to ( |22b| ) for /c = 
interpolating from say 6*1 = at a: — >• — oo to ^2 = 7'"/2 at a; — > +00 will thus give us its 
complete behavior: 

e = - arcsin(tanh(^)) + ^. (23) 

Thus, the near false vacuum state around the bubbles behaves as a rather effective insulating 
material as far as conduction of phase waves is concerned, and we should not expect any 
significant propagation of the phase until the bubbles run into one another. It may also be 
worth noting how the phase interpolates from one of its values to the other by means of a 
wall-like structure similar to walls created by the modulus of the field. The thickness of this 
phase wall, 6p^ seems to be closely related to the thickness of the bubble wall 6m, although 
this relation will obviously change as the phase difference between the bubbles change. 

Once the bubbles collide and their walls begin to merge, the phase of each bubble will 
start propagating into the other one at the speed of light. While this merging is taking 
place, its net effect will be to widen the phase wall, so that its final thickness will be at least 
twice the thickness of the bubble walls for relativistic bubbles. 

Figure 1 shows the results of simulations of this process. We have plotted the phase and 
the bubble walls for a 1-dimensional, 2 bubble configuration with an initial phase difference 
of 7r/2, for three different times of the evolution: when they are about to start colliding, at 
the middle of the merging process, and at the point of its completion respectively. We can 
clearly see how the phase wall smoothens out and thickens while this takes place. 



B. Phase propagation inside the bubbles 

Once the bubble walls have completed their merging, the phase is free to propagate in the 
resulting single cavity. Since, especially in the thin wall case, the modulus of the field inside 
the cavity remains essentially constant and equal to its true vacuum value, the equation of 
motion for the phase simplifies to a wave equation 

d^d^'e = 0. (24) 

Kibble and Vilenkin have studied the process of wave propagation for the Abelian 
gauge model without dissipation. In the case they studied the bubbles move essentially at 
the speed of light, and the problem has a high degree of symmetry: if the axes are chosen 
so that the bubbles nucleate along the z axis at say (0, 0, 0, ±R), the whole bubble collision 
process will be invariant under the 3-dimensional Lorentz group S0(l,2) in the {t,x,y) 
subspace. The bubble collision occurs then along the surface z = 0, t^ — — y'^ = , and 
for any point in that surface there will be a frame of reference in which that is the point of 
first contact. Thus, symmetry dictates that ^ be a function 

e = e{T,z), (25) 

where = t^ — — . In our case, the damped motion of the walls obviously breaks 
Lorentz invariance as far as the bubbles motion and collision is concerned. The points at 
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which the two bubbles connect to each other in the z = plane however will still move at 
a very high speed (in fact much higher than the terminal velocity of the walls), and, since 
the phase is not affected by the damping, S0(l,2) Lorentz invariance should still be a good 
approximation in our case. Inserting the ansatz (|2^ ) for 6 in (|2^) yields then the equation 
for the phase 

d^^e + -dr9 - d,e = 0. (26) 

T 

As initial setup we will assume that the bubbles have been nucleated at ±R in the z axis 
while the collision takes place at t = 0. For the initial conditions we will take the thickness 
of the phase wall to be negligible when compared to the other scales of the problem, that is, 

e\r=o = eoe{z) , 9.^1^=0 = 0, (27) 
where e is the step function. Equation (p6| ) can then be solved to get 

{00 for z > y/t'^ — — y'-^, 
^o^7P==;f for \z\ < Vt^-x^-y^ (28) 
—6*0 for — z > y/t'^ — — y'^. 

This approximation will obviously work better for bubbles that move with relatively high 
terminal velocities. Firstly of course because the Lorentz invariance approximation will be 
more adequate, but also because faster wall motion leads to smaller phase wall thicknesses, 
as seen at the end of the last subsection, so that (|27|) is a better approximation for the 
phase wall. It will also clearly work better in the central region of the bubble than close to 
its walls, where the interaction between these and the phase wave, as well as the detailed 
way in which the connecting circumference between the two bubbles expands, may be of 
importance. 

To see how well this approximation holds we can compare it with numerical simulations, 
which will however be performed in 2 + 1 dimensions. The solution to the 2 + 1 version of 
is easily found 

{^0 for z > y/f^ - x^, 
2^ arcsin(^7^^) for \z\< VW^^, (29) 
-^0 for - z > y/f^ - x^, 

and has general features similar to those of the 3 + 1 case. Figure 2a shows a couple of 
snapshots of the phase wave profile along the z axis at different times after the collision for 
thin wall bubbles with relatively low friction. As we can see, our results give a rather accurate 
picture of how the phase evolves in this case. Figure 2b shows the same situation for the case 
in which the bubbles move under high friction (i.e., low terminal velocity). As we can see, 
here the phase wall thickness at the the time at which the bubbles finish merging is rather 
large, and therefore the step function approximation in (|2^ ) fails to correctly represent the 
initial state of the phase. We should note however that the phase propagation seems to 
proceed along similar lines as before, so that if we were to solve ( ^6]) with a smoother ansatz 
for 6{t = 0) we should again reproduce the observed phase behavior. 
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C. Interaction between the phase wave and the bubble walls: phase bouncing 



The most salient feature of the problem of bubble collision in the damped regime will 
be the possibility of interaction between the phase waves that propagate inside the merged 
bubbles, and the bubble walls. Such a situation was of course never encountered in the case 
of walls moving in vacuum and reaching asymptotically the speed of light. 

We can roughly anticipate the outcome of this interaction on physical grounds: since the 
phase is massless inside the bubbles, but massive in the near false vacuum outside them, 
it should follow that only the contributing modes to the phase wave which have sufficient 
energy to acquire the required mass will be able to go through the wall at all, while all the 
others must bounce off it. Whether these modes will exist at all in the phase wave, and if 
they do in what proportion, will then determine how the phase wave will behave after the 
bouncing. 

We will first study the interaction between the phase wave and the bubble wall both 
analytically and numerically in the simpler 1 + 1 dimensional scenario, and then try to 
extrapolate our conclusions to the 2 + 1 case and confront them with more numerical simu- 
lations there. 

In 1 + 1 dimensions, our damped thin wall solution looks like 



Ptv 

2 



1 — tanh 



1 



(30) 



Boosting to a frame of reference that moves along with the wall and has its origin at its 
center leads then to 



P 



Ptv 



(31) 



We will be interested in the situation in which phase waves approach the wall from 
X —oo,t — oo and look at the asymptotic outcome for t +oo, assuming that the 
incoming phase waves carry a very small amount of energy as compared to that stored in the 
wall so that we can neglect any back reaction on it due to the collision. In these conditions 
we can then take the shape of the wall as being essentially fixed by (|3T|). The equation of 
motion for the phase simplifies to 



e -9" - 2^9' = 0, 
P 



(32) 



or, doing a Fourier transform in t, 



9" + 2^9' + iu^9 = 0. (33) 
P 

Using then (|3TD for p as advertised and performing the change of independent variable 



\ _|_ Q—\/Xr)x 



(34) 
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puts the phase equation in the form 



y{l-y) 9" + (1 - 4y) 9' + ^^.^ _ = 0, (35) 

where now the prime indicates derivatives with respect to y. Equation (|35|) is a form of 
Papperitz equation with regular singular points at 0, 1, +oo, (see Ref. [jT^) and consequently 
its solutions will be given by hypergeometric functions. 

To learn about the fate of an incoming wave we will look at the two independent solutions 
of (|35D around y —>■ 1 (i.e., x +oo) and select the one that represents an outgoing wave. 
In terms of x, these outgoing wave solutions are found to be {B being an arbitrary constant) 



with uj= Y 1 — ^^^/(Ar]^), and where for notational convenience we have used 6m = 
for the bubble wall thickness. (Note that the exponential growth of the solutions would 
apparently yield a diverging energy density, specially for the case \uj\6m < 1, more on this 
later) . 

Using now the connection formulas for the hypergeometric functions, we can express 9 
in terms of incoming and outgoing waves around x — >■ — oo. We get 



9.,.^^^-(-ir-'B 



Fie^"^^ + Tae"*^^ . (37) 



where ri,r2 are combinations of gamma functions. Demanding that the incoming wave 
has unit amplitude fixes the constant B and determines the reflection and transmission 
coefficients. We obtain as final expressions for 9 

~ e"^^ + ^e-^-^ (38a) 
J- 1 

[ (Fi)"^ e^/^- e^'^/^™ ,\uj\6m < 1 

It only rests now to write down the precise form of the gammas and analyze their behavior. 
We get 

_ r(l + 2iuj6m)T{2- uj -iuj6m)Ti-l- ^ -iuj6m) . s 

r — - - v^yj 

^1 r(l - 2iuj6m)T{2- UJ +iuj6m)T{-l- oj +iuj5ra) 

Thus, for as long as uj remains real (i.e., for as long as Ic^l^m < 1), numerator and 
denominator in (|39D will be complex conjugate of each other, and the refiection coefficient 
will be 



R 



r 



2 



1. (40) 
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When \uj\5m > 1, uj^ i u however, and R is 



Sinh^ ( TT ( Uj5rr, — OJ 



R = ) ) (41) 

sinh TT u5m+ ^ 



and we see that i? ^ for u;(5m ^ 1. Thus incidents waves with wavelength Ijio larger than 
the wall thickness will be completely reflected by the wall, whereas waves with wavelengths 
shorter than the wall thickness will be partially reflected and partially transmitted, the 
reflected part going to zero quickly as the wavelength decreases. Of course these conclusions 
will have to be somewhat modified to describe the process from the point of view of a static 
observer at the origin. For such a reference frame, a boost with velocity —Vter yields for the 
new frequencies lo' of the incident and reflected waves 

J,^=^Uj{\\Vter). (42) 

, 1 - Vter 



Uj'^^f = (1 - Vter) 



Co",-, 



where (3 = l/yl — f|ej.. Thus, for a static observer the condition uj6m < 1 is transformed 
into 

U'ij^ < (3 {1 + Vter) , (43) 

while the reflected wave will present a Doppler shift relative to the incoming wave due to it 
having bounced off a moving wall. 

Figures 3a, 3b show an incident wave train of wavelength roughly five times that of 
the wall (in the static frame) being completely reflected. The reflected wave has the same 
amplitude as the incident one, but its frequency has decreased due to the Doppler shift. 
Figures 4a, 4b then show how a wave train with wavelength about half the wall thickness is 
partially reflected. Both cases present exponential growth of the phase outside the bubble 
(remember that, modulo 27r, the maximum phase difference that one can have is vr). This 
seems to fit the general behavior (p6|). Note however that an exponential growth of the 
phase with x does not neccesarily lead to an equivalent growth of the energy stored by the 
phase gradient, since this goes like ■i]'^ {86/ dx)"^ and the modulus is decaying exponentially. 
Thus, in the expression for the case uj6m > 1 the energy in the phase gradient will tend 
to a constant for instance. If the behavior of the phase outside the bubble were to be like 
the one given by the utSm < 1 regime in the same equation however, the phase gradient 
energy would diverge exponentially with x. Any of these two behaviours however will end 
up in the breakdown of our approximation that the wall is unaffected by the phase wave 
propagation. Once the phase wave gets through to the false vacuum, the modulus of the field 
stops behaving like a transparent medium as we saw in section III.l. In those conditions, 
and since the modulus decays exponentially, no matter how small was the energy carried 
initially by the phase wave we will always get to a point in which it will be of the same 
order than the energy of the modulus. At that point, backreaction on the modulus is no 
longer negligible, and a proper analysis of the problem would require solving the two coupled 
equations. The situation is depicted in Figures 3c, 4c, where a blow up shows how the wall 
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develops what appear as small bumps in its decaying region. Figures 3d, 4d then, where 
the energy of the modulus and of the phase are plotted, show how these bumps appear at 
the points where the energy of the phase is getting larger than that of the modulus. The 
result of this backreaction seems to be that the modulus absorbs the excess energy in the 
phase, whereas the phase will still propagate forward in some exponential way. We have not 
carried out a detailed analysis of this part of the problem however, since it does not directly 
affect the behavior of the phase inside the bubbles in any relevant way. 

To sum up, we can expect that waves with wavelengths (roughly) larger than the bubble 
thickness will be totally reflected by the wall, the rest being partially transmitted and 
partially reflected by it. It is then easy to see what will happen to the phase wave in a (1 + 1 
dimensional) two bubble collision: as we saw in section III.l, the phase wall thickness at 
the end of the collision will be at least of the order of twice the bubble wall thickness, and 
more likely larger than that. It seems clear then that basically all of the Fourier components 
of the phase wave will have wavelengths that will fall into the total reflection regime, and 
thus the whole phase wave itself will simply be reflected by the wall, while some form of 
exponential behaviour propagates to the outside of the bubbles. Let us imagine for the sake 
of clarity that, at collision time, bubble 1 had a phase 6i and bubble two 62 > 61, with 
02 — di = A^. If the shape of the phase wall at the completion of the collision was f{x) then, 
after it, we will have phase waves with shape f{x — t)/2, f{x + 1)/2 propagating into each 
bubble, carrying a phase difference —A6/2 into bubble 2 and A6/2 into bubble 1. After 
these waves have bounced off the bubble walls and propagated back into the interior again, 
the phase of each bubble will however be, for bubble 2, 62 — 2A6/2 = 61, and for bubble 1, 
61 + 2A6/2 = 62- The phases of the bubbles will thus have switched. The whole process is 
depicted in Figure 5, where the refered sequence has been plotted from a simulation. In Fig. 
5a, the walls of the two bubbles are just about to finish their merging (continous line), and 
the shape of the phase wall at that time is shown (dashed line). The bubble to the right 
plays the role of bubble 2 above, having 62 > 61. The following pictures show how the two 
phase waves propagate into the bubbles (Fig. 5b), and back after bouncing (Fig. 5c). As 
expected, after reflection each phase wave still carries a phase of ±A6/2, for bubbles 1 and 
2 respectively. Finally, in Fig. 5d the two phase waves meet again. The phase polarity of 
the system has been completely inverted. 

If this mechanism had no energy losses, the two bubbles would then behave as a sort 
of resonator in which the phase of each bubble would continuously oscillate between 61 and 
6*2, their respective polarities always switched. Since the bubble walls are moving though, 
the phase waves present the expected Doppler shift after bouncing off the walls as was seen 
before (which can also be clearly seen in Figs. 5c, d). For non-relativistic terminal velocities 
(i.e., for Vter ^ 1), the magnitude of this shift is from (^2]) (we drop the primes in the 
notation for the frequencies here) 



or what is the same, after n reflections, the resulting phase wave will have a frequency 



The oscillation process will effectively die off when ujnref ~ (2"^^) ^, 'R- being the (approxi- 
mate) radius of the single true vacuum cavity after n reflections have taken place. If uj~^ is 



UJref ~ CJi„(l - 2Vter) 



(44) 



IVnref ^ UJin{l - 2VterT ■ 



(45) 
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very small when compared to TZ, and neglecting higher orders in vter, ^nref becomes then of 
the order of the radius after 

n ~ ^— (46) 

oscillations. On the other hand, it is not difficult to see that the time needed to complete 
each oscillation grows like 

t^ = ARj-^-X, (47) 

V 1 - Vter J 

where Rq is the radius of the bubbles at collision time (keeping in mind that it will take 
only half tito reach the walls for the first time though). Thus, the relaxation time until the 
oscillations die off will be 

9 F? n~l/2vter , 1 ^ n 



For a terminal velocity around 0.1, (^Uj) gives us five oscillations until equilibration, for a 
total oscillating time Tosc ~ 24i?o from (^8]). We see therefore that the relaxation time for 
the system can in general be quite significant. 



IV. THREE BUBBLE COLLISIONS: VORTEX FORMATION 

We wish now to try to generalize the results found for phase wave bouncing in 1+1 
dimensions to 2+1 dimensions, since our final concern is to find whether these effects can 
affect the formation of vortices. Finding an exact solution for the phase propagation and 
interaction in 2+1 dimensions is an extremely involved problem, but the results in one spatial 
dimensions can certainly be used as a guide, and numerical simulations can provide the rest. 

Note first that wave propagation in two or three spatial dimensions differs qualitatively 
from the one dimensional case. Whereas in one dimension we have two phase waves that 
propagate inside the bubbles carrying half of the phase difference each, in two or more 
dimensions we have two wavefronts with amplitudes that decrease in time and a region that 
continuously interpolates between them (see the previous section). In Figure 6 we have 
plotted magnitude and phase contours in a two-bubble collision, as well as the phase at 
each point represented by an arrow. We can see how the phase interpolates between the 
two bubbles at merging (6a), and then starts to propagate. During the initial stages of 
propagation, the points of contact of the two bubbles move in the direction perpendicular to 
the z axis at superluminal speed, so that the phase waves cannot reach them. Propagation 
is similar in these first stages to the undamped case: the region of interpolated phase is just 
a semicircle (6b). Soon however the contact points will slow down to asymptotically reach 
the terminal velocity. When the phase waves reach the walls, the shape of the interpolating 
region is affected by the interaction, as can be seen in 6c and d. We come then to the first 
consequence of the damped motion of the walls. The fact that the phase waves are able to 
catch up with the bubble walls means that, at this point, a region of smoothly interpolated 
(i.e., nearly homogeneous) phase extends up to the boundaries of the system. At this stage 
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then it is easy for a third bubble of the right size and position to colhde with the merged 
two-bubble system inside the region where the phase is nearly homogeneous. Hence, no 
vortex will be formed. The situation is depicted in Fig. 7, where a third bubble collides 
with the two bubble system while the processes seen in Figure 6 are taking place. In 7a 
and b, the third bubble collides shortly after the phase starts interpolating, so that a vortex 
is created. In 7c and d however the nucleation of the third bubble occurs a little later. 
When it finally reaches the two bubble system, the phase has already interpolated beyond 
the collision points and no vortex is formed. Notice that the initial phases of the bubbles 
in these figure are not exactly 7r/6, 57r/6 and 37r/2, but the two initial bubbles have slightly 
smaller values of the phases, so that the resulting interpolated phase after merging is slightly 
smaller than 7r/2. This was done in order to avoid having the third bubble colliding with 
a region of exactly opposite phase. As was shown in Ref. [0, such a collision produces a 
long-lived "domain wall" state between the bubbles that delays the merging significantly. 
We will come back to this point after we discuss the further evolution of the system below, 
for the time being we just remark that the phase distribution is responsible for the slight 
asymmetry in the vortex's position. 

The situation depicted in Fig. 7 could represent a rather efficient mechanism for the 
suppression of vortex formation, if it were not for the interaction with the walls. From the 
1-1-1 case, one expects the phase wave to bounce back, thus spoiling its homogenization. 
How does the bouncing look like in 2-1-1 dimension? In two spatial dimensions, at any 
point that the phase wave meets the wall, its propagation vector will have one component 
in the direction normal to the wall and another one tangential to it. Only the component 
in the normal direction will see the wall and bounce off it however, while the other one will 
continue to propagate freely. Thus, although in general the interaction between the wall and 
the phase wave will be a complicated superposition of these two processes, we can expect 
that after some time, in the region close the walls the phase will predominantly propagate 
tangentially to them, the rest of it having bounced towards the center. The effects of this 
bouncing can already be seen in Figures 6c and d: it is the reflection of the normal modes 
that causes the change of shape in the wavefront. 

After this we only have to wait for the central region of the phase, propagating along the 
z axis, to get to the end of the bubble and bounce off the wall. Fig. 8 is a continuation of the 
same simulation started in Fig. 6. At the stage depicted by 8a, the phase wave has reached 
the end of the bubble. Therefore, at that point virtually all the phase propagating normal 
to the wall has bounced towards the center, and tangential propagation dominates close to 
the wall. In Figure 8b we see the subsequent evolution of the system. While the central 
region of the phase, propagating along the z axis, collides head on and bounces off the end 
of the bubble as expected, the tangential components cross each other at that region and 
start propagating in the opposite direction, back towards the center of the collision. Thus, 
the combination of these two phenomena brings about an inversion of the phase similar to 
that found in 1-|-1 dimensions. Finally, a third bubble collides with the system in 8c, after 
the bouncing has taken place. Because of phase inversion, we get an antivortex formed in 
8d. 

Note that if we were to wait for the collision to happen until yet another bounce had 
taken place (for a total of two bounces), the resulting defect would again be a vortex. 
Three bounces and an antivortex again, and so on. To get a total suppression of the vortex 
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formation probability, therefore, one would have to delay the collision of the third bubble 
until the system has relaxed. As we have seen, the 1+1- dimensional model predicts a large 
value for the relaxation time. Simulations show that for Vter ~ 0.1, vortices can still be 
formed in collisions occurring after the merged system has gone through five oscillations, 
our predicted time for equilibration. So the two-bubble system will have to remain in 
"isolation" for a very long time to be able to acquire an homogeneous distribution of phase, 
a very unlikely situation. Our results seem to indicate then that an analysis of the vortex 
formation probability will have to take into account multiple bubble collisions, such as for 
example the one carried out in Ref. 

The situation is of course further complicated if we consider different initial positions 
for the colliding bubbles. We have taken them to be situated so that the latest bubble is 
centered in the plane of collision of the other two, but it is easy to see that any shift in 
position will affect the formation of vortices in a non trivial way, as it does any change in 
the nucleation time. We conclude with some more words about the case depicted in Fig. 7c 
and d. Had we set the three phases equally distributed, the delay caused in the merging due 
to the formation of a metastable wall between the region of opposite phases will have given 
enough time for the bouncing to take place, and an antivortex to be formed. Simulations 
show that for the positions of the bubbles considered here, a defect is always formed in 
the particular (and less likely) case of equally distributed initial phases, unless the phase 
oscillations have completely relaxed. 



V. CONCLUSIONS 

We have analyzed bubble motion and interactions in a plasma, simulating its damping 
effect by a friction term in the equations of motion. We have found that there exists an exact 
solution for the damped motion of the bubble in the thin-wall regime, representing a bubble 
propagating with a terminal speed roughly equal to the inverse of the friction coefficient. 

An analytical study of the collision of two bubbles, and the subsequent process of phase 
interpolation, was then performed. First, the shape of the "phase wall" that interpolates 
between the bubbles was found, and it was shown that the vacuum state outside the bubbles 
is a very efficient "insulator", so that one does not expect the phase to propagate in the 
false vacuum. The initial thickness of the phase wall was shown to be at least twice that of 
the bubble's walls. 

Once the initial set up for the collision was thus determined, the equations for the phase 
propagation were solved, and the interaction with the walls determined, using the previously 
found solution for the wall's motion. It was shown that phase waves generated in a collision 
have wavelengths such that they are always reflected by the walls. The result is that an 
oscillating state is formed, in which the bubble's walls act almost as a resonant cavity. Since 
the true vacuum cavity expands, though, these oscillations are damped and eventually die 
off. For terminal velocities typical of those expected in an electroweak phase transition, the 
relaxation time is estimated to be around 2AR, with R the radius of the bubbles at collision. 

Using numerical simulations, these results were shown to hold in 2+1 dimensions. The 
reflection of the phase wave goes along similar lines than the one-dimensional case. The 
formation of vortices was shown to be affected by the damping motion, basically due to 
two facts: flrst, that the interpolated phase region can reach the boundaries of the system. 
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where a reflection of the phase waves ofi^ the bubble walls will take place; and second that 
the long-lived oscillatory state thus produced delays significantly the homogenization of the 
phase. The main effect of this oscillatory state will have on vortex formation processes will 
be that it will become possible, for the same set of three bubbles, to form a vortex, and 
antivortex or no defect at all depending on the precise timing of the last collision. 
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FIGURES 



FIG. 1. Modulus of the field (a) and phase (b) of a system of two bubbles, along the axis that 
joins their centers, for three moments of the collision process. Continuos line is for t = 16, dotted 
for t = 19, dashed for t = 20 (where radial and time coordinates are given in units of the field's 
mass, and the modulus is normalized with the symmetry-breaking scale). Bubbles are nucleated 
at t = 0, and 7 = 3. Here and in the following graphs, e = 0.8. 

FIG. 2. Phase propagation inside the bubbles. In (a), the friction coefficient 7 is 3, in (b) is 10. 
The analytical solution is shown for different times with a continuos line. Results of the simulation 
of the phase propagation are shown with triangles (for t = tc + 1, tc= collision time ) and circles 
(for t = tc + 7). 

FIG. 3. Interaction of a phase wave with the bubble wall, for a phase wavelength approximately 
equal to five times the width of the wall. Arrows indicate the direction of propagation of bubble 
and phase. In (a), (b) and (c), the field's phase (dashed line) and magnitude (continuos line) are 
plotted for different times as indicated. In (d), the energy in modulus (continuos line) and phase 
(dashed line) are compared. 

FIG. 4. Same as Fig. 3, where now the phase wave has a wavelength of approximately half the 
wall's width. 

FIG. 5. Phase propagation and bouncing. The phase (dashed line) and the bubble walls (con- 
tinous) are plotted for {a)t = 15, (h)t = 25, (c)t = 50 and {d)t = 60, with 7 = 10. 

FIG. 6. Collision of two bubbles with phases 7r/6 and 57r/6, and subsequent phase interpolation. 
Continuos lines axe contours of equal modulus of the field, dotted lines are phase contours. The 
phase of the field is also represented by arrows. In this and the following graphs, 7 = 10. The axes 
are arbitrary lattice coordinates. 

FIG. 7. Collision of three bubbles with phases {tt/6 - A), 57r/6 - A) and 37r/2, with A = 0.01. 
(a) and (b) are plots of two moments of the collision of two bubbles nucleated at i = with a third 
one nucleated at t = 8. (c) and (d) correspond to a similar simulation, where now the third bubble 
is nucleated aX t = 18, in the same position as the previous one. Only field magnitude contours 
are represented. 

FIG. 8. Later evolution of the system of Figure 6. A third bubble is nucleated at i = 48. Phase 
contours are omitted in the last frame for clarity. 
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